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. Abstract 



We generalize the Lewis- Riesenfeld technique of solving the time-dependent Schrodinger equation to cases 
where the invariant has continuous eigenvalues. An explicit formula for a generalized Lewis- Riesenfeld phase 
is derived in terms of the eigenstates of the invariant. As an illustration the generalized phase is calculated 
for a particle in a time-dependent linear potential. 
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The study of time dependent quantum systems has attracted considerable interest in the litterature. The 
origin of this development was no doubt the discovery of an exact invariant by Lewis and Riesenfeld [1 . 
The work of Lewis and Riesenfeld and others assumes that the eigenvalue spectrum for the invariant I is 
discrete. Let us recall that the general method to introduce the Lewis and Riesenfeld theory, valid whatever the 
time dependence of the parameters, considers invariant operators. For a system specified by a time-dependent 



Hamiltonian H ( X (t) ) , and a corresponding evolution operator U (t), an invariant is an operator I(t) such that 



> 
On 
OO 

I(t) = U(t)I(0)U- 1 {t). (2) 

It possesses a remarkable property that any eigenstate of 1(0) evolves into an eigenstate of I(t). Then, if the 
set of reference eigenstates {| <j> n (t)}} for the operator I(t) are continuous with respect to t (all eigensates are 
associated with the same time-independent eigenvalue e n ), the corresponding global phases 6 n {t) are defined 
by the relation associated to the wave functions \ip n (t)): 



K(t)) = [/(t)|0„(O)) =e ie ^\cj> n {t)). (3) 
^ . It follows from the Schrodinger equation for \tp n (t)) 



that n {t) satisfies the relation 



^ \1>n{t)) = H(t)i> n (t) K(f)> , (4) 

n^9 n (t) = (tf, n (t)\ih^-H\4> n {t)). (5) 

One way to describe the exact quantum evolution of Lewis and Riesenfeld is to introduce the concept of 
elementary projectors on an eigenstate | 4> n (t)} of the invariant operator I(t) 

P n (t) =| 4> n {t))(<t> n (t) I • (6) 
It is easy to verify that each projector P n (t) is therefore a constant of the motion i.e., P n (t) = U (t) P n (0) U + (t) . 
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We can state that the exact evolution described by equations (3) and (5) can be formally written in the 
following form 



Vt: U(t)P n {Q) = P n (t)U(t). (7) 

Notice that if, initially, the system is in the eigenstate \<j) n (0)) so that / (0) \<fr n (0)) = e n \(f> n (0)), then 
P n (0) |0 B (0)} = \<t> n (0)> and © gives 

ff(i)|0 n (O)) = P n (t)tf(t)|^ n (O)>. (8) 

In general the spectrum of / possess both discret and continuous eigenvalues. The Lewis and Riesenfeld 
theory in a continuous spectrum was raised for the first time by Hartley and Ray [2] where they extend the 
Lewis and Riesenfeld theory for a general Ermakov system to the continuous spectra like an anstaz and looks at 
the eigenfunctions in a continuous spectrum \cf>(k;t)) of the invariant operator I(t) and the solution \%l)(k,t))oi 
the Schrodinger equation in the form e 1 " 9 **^ \<p (fc; t)) .The limitation of the Hartley-Ray approach [5] is that, in 
general, there is no explicit formula of Lewis and Riesenfeld phase. Later Gao et al [3] calculate the Schrodinger 
solutions for a continuous spectrum using the path integral technique. In some papers [3] [3] the overall phase 
factor 9k(t) , interpreted in the spirit of the original investigation of Lewis and Riesenfeld, is actually obtained 
through the relation h-^9k(t) = (c/)(k;t) \ ih-^ — H | 0(fc;i)).This procedure clearly parallel to that for the 
discrete spectrum is not founded. The reason is that the calculation of the expectation value (ih-§[ — H) with 
respect to the eigenfunctions \<p{k]t)) analyzed trough examples in [4| lead to errenous results. 

In the case of continuous spectrum we cannot numerate eigenvalues and eigenfunctions, they are characterised 
by the value of the physical quantity in the corresponding state. Althoug the eigenfunctions \<j>{k\t))oi the 
invariant operators with continuous spectra cannot be normalised in the usual manner as is done for the functions 
of discret spectra, one can construct with the \<f){k]t)) new quantities - theWeyl's eigendifferentials (wave 
packets )- [5] [6] which possess the properties of the eigenfunction of discrete spectrum. The eigendifferentials 
are defined by the equation 

k+Sk 

\8<j> (fc; t)) = J \<p(k';t))dk'. (9) 

k 

They divide up the continuous spectrum of the eigenvalues into finite but sufficiently small discrete regions of 
size Sk . 

The eigendifferential ^ is a special wave packet which has only a finite extension in space; hence, it vanishes 
at infinity and therefore can be seen in analogy to bound states. Furthermore, because the Sip have finite spatial 
extension, they can be normalized. Then in the limit Sk — * 0, a meaningful normalization of the function ip 
themselves follows: the normalization on S functions. 

For Sk, a small connected range of value of the parameter fc (this corresponds to a group of "neighboring" 
states), the operator 

k+Sk 

SP(k;t) = J \<j}(k';t))(<j)(k';t)\dk' (10) 

k 

represents the projector (the differential projection operator [5 6 ) onto those states contained in the interval 
and characterized by the values of the parameter fc within the range of values Sk. The action of SP (fc; t) on a 
wavefunction \ip (t)) is defined by 

fc+<5fc 

SP{k;t)\tp(t)) = J C(k';t)\cl)(k';t))dk'. (11) 
fc 

The application of the differential projection operator SP (fc; t) causes thus the projection of the wavefunction 
onto the domain of states <fi (fc; t) which is characterized by fc values within the Sk interval . 

In this letter, we present a straightforward, yet rigorous, proof of the exact quantum evolution for systems 
whose invariant has a completely continuous spectrum supposed to be non-degenerated for reasons of simplicity. 
Example of a particle in a time-dependent linear potential is worked out for illustration. The case of both disrete 
and continuous eigenvalues can be obtained by superposition. Before proceeding further, we give the statement 
of the exact quantum evolution. 
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Given a physical system with a time-dependent Hamiltonian H (t) , it is possible to build an invariant operator 
I(t) verifying |T]), such that its eigenvalues Sk are purely continuous and constants 



J(t)|0(fe;t)>= Efc |0(fc;*)>. (12) 

Let us call U (t) the evolution operator associated to the time-dependent Hamiltonian H (t). The evolution of 
the system obeys to the Schrodinger equation 

ih^m))=B{t)m)), (is) 

under these conditions it is possible to state the exact quantum evolution: 

" If the quantum system with time- dependent Hamiltonian whose invariant has a completely continuous spec- 
trum supposed non-degenerated is initially in an eigenstate \4>{k, 0)) of I (0) then the state of the system at 
any time t will remain in the subspace generated by the eigenstates \(f>(k;t)) of I(t) pertaining to the interval 
[k,k + Sk} , \ 

In others words, the exact evolution can be formally written in terms of the evolution operator as 

Vk,Vfc: U(t)5P(k;Q) = 5P(k;t)U(t). (14) 

The proof of this last equation is straightforward. One has to take the partial derivative of SP(k;t) — 

k+Sk 

J \<p (k'; t)) {4> [k'\ t)| dk' with respect to the time t, taking into acount Eqs. (1) and (12), one will get then: 

k 

^2^ + ^[SP(k;t),H]=0, (15) 
at in 

which has a formal solution SP (k; t) = U (t) 5P (k; 0) U+ (t) . 

Notice that if, initially, the system is in the state \ip(k, 0)) = \<p(k,0)), then (fHf| implies that 

SP (k' ; t) U (t) |i/> (fc, 0)} = U (t) 5P (k'; 0) \ip(k, 0)) , (16) 

expanding an arbitrary state vector (t)) on the basis of the instantaneous quasi-eigenfunction of I(t) and 
using (11) we obtain 

fc+5fc 

|V (*,*)>= J C{k'-t)\<j){k'-t))dk'- We[k,k + Sk], (17) 

'k 

we see that the state \ip (k, t)) belongs to the subspace generated by the states \4> (k; t)} pertaining to the interval 
[k, k + Sk}. 

Inserting (flTf in the Schrodinger equation (|T5)) . lead to 

k+5k k+Sk k+5k 

ih^C (k';t)\(/)(k';t))dk' + J ihC (k 1 '; t) ^ \(f> (k 1 '; t)) dk' = J C (k';t) H (k';t)\ip (k';t)} dk' . (18) 

k k k 

We multiply Eq. (TTH| by the bra of the eigendifferential Q introduced earlier 

7+<57 

(6<f>( T ,t)\= J (4>(v,t)\dr) (19) 

7 

this yields 

k+Sk k+Sk 

ih—C(k';t)\(t)(k , ;t))dk'= / C(k';t){S(f){r,t)\H{t)-ih—\(t)(k';t))dk' (20) 

k k 
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Since k can sweep all the possible values and the intervals 6k should be small (5k — > 0) , the equality (|20|) between 
integrals implies the equality between integrands, hence 



ih^-C(k';t) = C{k'-t) 



This equation is easily integrated and gives: 



C(k';t) =5{k' - k)exp 



hence 



<(k;t)\H(t)-ih^- t \c/>(k';t)) 



: k' e[k,k + Sk] . 



(21) 



.(*;0|itf(0-^ltf(* , ;* , )>)*' 



; k' E [k,k + Sk], (22) 



(23) 



where 9 k (t) is the global phase given by 



h-U)= I {5^{k ] if)\i^j-'^H{i;)\<i } {k;t'))d£. 



(24) 



This explicit formula of the phase could not be made in the Hartley-Ray approach [2] as mentioned earlier. 

As the interval [k, k + Sk] is located inside of the interval [—00, +00], we can write the generalized phase in 
the following practical form 



t +00 



*(*) 



. (k'; t') I % A -iff (A; t'))dt'dk' 



(25) 



-00 



which embodies the central result of this paper. 

To illustrate this theory, let us calculate this phase for a particle moving in a time dependent linear potential 



2m 

where f(t) is a time-dependent function. We look for the invariant of the form 

I(t) = a (t) p 2 + b(t)p + c{t)x + d (t) . 
The invariant equation (1) is satisfied if the time-dependent coefficients are such that 

a = ao, 



(26) 
(27) 

(28) 
(29) 
(30) 

(31) 

where ao, 60 1 Co and do are arbitrary real constants. We can choose ao = 1 and do = without loss of generalities. 
The eigenstates of I(t) corresponding to time-independent eigenvalues k are the solutions of the equation 



b = 2a [ f dt' -c [ —dt' + b 
Jo Jo m 



c = c , 



d = 2a / / / / dt'dt" -co f 



dt'dt" +b f dt' + d , 



d 



% — — ^ — ihb— h cqx + d 



cf> k (x,t) = k(p k (x,t) . 



(32) 



dx 2 dx 

The key point of our analysis is to perform the time-dependent unitary transformation such that 

$ k (x)=S(t)4> k (x,t), (33) 
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where a time-dependent unitary operator H(t) is given by 



1 ; 7^(t '') /f ) x " x,) (4- r i 



(34) 



It can be easily shown that, under this transformation, the coordinate and momentum operators change ac- 
cording to 



S(t)xS(t)+ = x + 1 (j- - i ) . 
— > H(i)pH(i)+=p-^. 



(35) 
(36) 



An important property of the transformation exp y— i-^ ^ — dj pj , the action of which on a wave function 
in the x representation reads 

1 (b 2 



exp 



- l k (r- rf H /, '- /! = / 



x I — — d ] ,t 

c V 4 



(37) 



Hence, the operator / changes into time-independent operator Iq = S / S + = p 2 + cox. In other words, the 
eigenvalue equation (32) for the transformed invariant operator can be simply represented in the form of Airy 
equation 

d 2 



dZ 2 



*k(Z) =0, 



where we have introduced a new variable Z related to x through the relation Z = (-p-) 3 (a 
of the fundamental one dimensional ordinary second-order differential equation (38) is well-known 



(38) 

r-).The solution 



* f ^ 1 ( 1 



2ir \c h 4 



k 

x 

Co. 



(39) 



+ 00 



Ai (x) — J e lXx e l 3 d\ being the integral representation of the Airy function. It is easy to verify that 

—00 

l$fe> = 6(k-k'). 
Reversing the procedure above, we can obtain 



4> k (x,t) 



1 / 1 



2tt \c h 4 



exp 



2h 



Ai 



£o\ 3 

h 2 . 



1 (b 2 , 

a; — + k — d 

co V 4 



(40) 



There remains the problem of finding the phases Ok(t) which satisfy (25). Carrying out the unitary trans- 
formation S the right -hand side of Eq. (25) becomes 



iH Wt -H(t) 



k / - 



1 

2m 

1 
2m 



$k> 



S(k-k'), 



(41) 



where we have used Jo l^fc) = k |3>fc) and the <5 Dirac normalisation of Note that if k = k' the latter Eq. 
(41) leads to an unfounded result (equal to infinity) as was mentioned earlier. Hence the phase (25) is found to 
be 



Z 

0k(t)-0 k{ 0) = -±J 



k+ b l-d 



dt'. 



(42) 
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